0-7T transition in SFS junctions with strongly spin-dependent scattering 
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We develop theory of the critical current in a superconductor - ferromagnetic alloy - superconduc- 
tor trilayers, which takes into account strong spin dependence of electron scattering of compositional 
disorder in a diluted ferromagnetic alloy. We show that in such a system the critical current os- 
cillations as the function of the thickness of the ferromagnetic layer, with the period of vp/2I, vf 
and I being the Fermi velocity and exchange splitting, respectively, decay exponentially with the 
characteristic length of the order of the mean free path. 
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The recent observation-^^ of Josephson junctions 
with negative coupling, 5,6 also known as ir junctions, has 
attracted a lot of attention to hybrid superconductor - 
ferromagnet - superconductor (SFS) structures. In con- 
trast to conventional Josephson junctions, such as su- 
perconductor - normal metal systems, where the ground 
state corresponds to the superconducting phase differ- 
ence if of zero, the phase difference in a SFS trilayer can 
take both tp = and ir values, depending on the thickness 
of the ferromagnetic layer. Both and tt states in SFS 
trilayers have been deduced from the measurements of 
the density of states^ and the critical current as a func- 
tion of magnetic flux and temperature.—^ 7 -^ In particu- 
lar, the critical current exhibits oscillations superimposed 
on the exponential decay as a function of the thickness 
of the ferromagnetic layer^ The decay length £<j and 
the period £ Q of these oscillations have been measured, 
providing comparable yet unequal experimental values of 
these two parameters. 

In ballistic SFS structures, the critical current is ex- 
pected to oscillate with the period £ = vp/2I. For a 
ferromagnetic metal with a strong exchange splitting /, 
fluctuations of the width of the ferromagnetic layer sup- 
press the appearance of the proximity effect, despite the 
fact that in ballistic structures Cooper pairs decay with 
the distance according to a power law rather than ex- 
ponentially. Moreover, it has been showni - that, when 
the electron motion in a ferromagnetic film with large 
I is diffusive, the randomisation of the oscillation phase 
over paths of different lengths leads to the expenential 
suppression of proximity at the length of the mean free 
path: ~ I for the case of /r > 1 (where r is the 
electron mean free path). To enhance the proximity ef- 
fect in a SFS multilayers, one may want to use weakly 
ferromagnetic alloys, where the exchange field / is re- 
duced by diluting the magnetic component. The analy- 
sis of diluted systems with It <C 1 based upon modelling 
disorder in SFS junctions as spin-independent impurities 
has shown that the decay length may be expected 6 -^ 
to extend beyond the mean free path range, such that 
£i ~ Co = y/D/l, where D = v 2 F r/3. 

In this paper, we show that a possibility to prolong 
the extent of the superconducting proximity effect in SFS 



structures by making them of diluted magnetic alloys is 
strongly limited. Following theory of suppression of su- 
perconductivity by magnetic impurities^, earlier theo- 
rie s 13 ' 14 took into account the effect of magnetic disorder 
by including in the Usadel equation a weak Cooper pair 
relaxation described by a phenomenological spin relax- 
ation rate r~ ■ Keeping in mind that even in a weak 
ferromagnet electron spin flip is an inelastic process and 
should be accompanied by the excitation of a magnon, we 
attribute the pair breaking in a ferromagnetic alloy to a 
giant magnetoresistance (GMR) type effect. As noticed 
in earlier GMR studie s 15 ' 16 , a feature of ferromagnetic al- 
loys is that elastic electron scattering in them is strongly 
spin dependent. Indeed, one scattering event off strongly 
spin-dependent disorder, seen differently by spin-up and 
spin-down electrons, is enough to break a singlet Cooper 
pair. In such a case, the decay length of a Cooper pair is 
of the order of the mean free path, £<j ~ I. Since, in this 
case, the use of Usadel equations adopted in the previ- 
ous studies of disordered SFS junctions 6 -^ does not hold, 
here we employ a nonlocal approach based on solution of 
Eilcnberger equation 6 - ' 10 ' 13 ! 17 to describe 0-7r Josephson 
oscillations as a function of the thickness of the diluted 
ferromagnetic alloy layer. 

To describe a dilute ferromagnetic alloy, we use the 
following Hamiltonian (a 2 x 2 matrix in the spin space) , 
adopted^ 5 - in GMR theory, 

H = p 2 /2m + V(r) + erJ(r) , (1) 

where V and J describe magnetic atoms embedded into 
a normal metal, and er is the vector of Pauli matrices. 
The average (J) = e z I determines the exchange splitting 
for conduction band electrons, and (V) = 0. Since ev- 
ery magnetic atom produces both scalar V and exchange 
J potentials, we use the following correlation functions 
for magnetic and nonmagnetic disorder, (V(r)V(r')) = 
{2^T V r 1 5{r - r'), (J a (r)J (r')} = {2mvrj)-H a pS{T - 
r'), and (V(r)J a (r')) = {2niyT mix )- 1 S az S(r - r'). 

The starting point for quantitative description is Eilcn- 
berger equation for the retarded component of the semi- 
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classical Green's function, 

v F nd r g + [-iuJT 3 + ila z + it, g] _ = 0, (2) 
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where (/ + ) Q/3 (r, t; n, w) = -[/ a( g(r,t; -n, -w)]*, r 3 acts 
in the Nambu space, n = p/p, the self-energy has the 
form 
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and (g) = J gd 2 n/4:TT is the Green's function averaged 
over momentum direction. For the weak proximity effect, 
Eq. ^ can be linearized around the zero-order Green's 
function go = t 3 . Performing the expansion up to first 
order, we obtain 

v F nd r f — 2iuif + 2ila z f+ 

(ry 1 +rj 1 )f-(ry 1 -rj 1 )(f) = 0. (4) 

The linearization of the Eilenberger equation and sub- 
sequent analysis are based upon the assumption of weak 
coupling between superconductors and the ferromagnet, 
which is realized, for instance, if these are separated by 
an opaque barrier with the low transparency 9 <C 1. 
The appropriate boundary conditions have been derived 
by Zaitsev^ 
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9s = 9f =9 > 

g a {(i - G)(l - (g a ?\ + Q{g s -) 2 } = &9 S -9% 

where g/ a — (gi(n z )±gi(—n z ))/2 (i — S and i — F for a 
superconductor and a ferromagnetic alloy, respectively), 
n z > 0, where n z is the projection of n = p/p onto the 
direction normal to the SF interface, and g± = (g^ ± 
g F )/2. In the case of low transparency <C 1, we find 
that in the first order in 0, 

where g)p and g^P are the Green's functions in the two 
materials when those are detached (0 = 0). Together 
with Eq. (U|) this gives us closed set of equations. 

It is convenient to represent the semiclassical Green's 
function /(n, z) as a combination of two functions of 
a positive argument, n z > 0: fi(n z ,z) = f(n z ,z) and 
f2(n z , z) = f(—n z , z). In this representation the bound- 
ary conditions take the form 

fi(n z ,0) - f2(n z ,0) = a L , 
fi(n z ,d F ) - f2(n z ,d F ) = -or, 



o-l/r = -©A cxp(i<t) L/R )/ VA 2 -u 2 , (5) 

where dp is the thickness of the ferromagnetic layer. The 
equations for /i and /2 take the form 

n z d z fi(n z ,z) + \fi{n z , z) - a(f(z)) = 0, 

-n z d z f2(n z , z) + Xf2(n z , z) - a(f(z)) = 0, 
where n z > and a = (tj — Tv) / {t~v + tj), —1 < a < 1. 
Also, 



is a 2 x 2 matrix acting on the 2x2 matrix /, and r 1 = 
Ty 1 + Tj 1 . The averaged Green's function equals 

(/(«))=/ dn z (h(n z ,z) + f 2 (n z ,z)). (8) 
Jo 

In the case of a thick ferromagnetic layer, such that 
e -d F /i ^ where 2 = v f t is the mean free path, one 
can write down the formal solution of Eqs. (|6|) as 

h(n z ,z) = a L e- Xz ^ 1 + j f" e x ^ '-'V"' 1 (f(z')) — 

I Jo n z 

e -A(*+»')/n.i^( z /)^ j ( 9a ) 
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e \(z'+z-2d F )/nJ^( z ,^^_^ (g b ) 
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The subsequent algebra includes adding and averaging 
Eqs. which leads to the integral equation for (/). 
Having presented (f(z)) as the sum, 



(f{z)) = a L h{z) + a R h(d F - z) 



(10) 



we find that the (matrix) function h(z) satisfies Fredholm 
equation of the second type, 

2/1(2) = K(Xz) + (a/l) J F (g(X\z - z'\) + 

G(\(z + z')) + G(A(2rf F - z - z')))h(z')dz', (11) 

where G{z) = J* n^e-^^dn^ K(z) = e - z / nJ dn z . 

Up to this point, we could still reduce our equations to 
Usadel equations provided the diffusion approximation 
holds, (1 — a)||ImA|| <C 1. In the rest of the paper, 
we work outside this regime and consider the ballistic 
situation. For a — 0, the exact solution of Eq. (fTTj) 
is h(z) = K(Xz)/2. Generalizing, we find that in the 
ballistic case the solution is determined by behavior of 
functions K and G which at z ^> I are K(z) « G(z) w 
e~ z l l l/ z. Assuming that solution falls off exponentially 
as e _Az /', one can see that in Eq. (fTTj) the last term in the 
integral can be neglected everywhere except for a small 
region near the boundary, z = d F - This enables us to 
split the solution of Eq. (|11[) into two parts, 

h{z) = e- Xz/l h L {z) + e- xt - 2dF - z ^ l h R (d F - z). (12) 

The first term is relevant everywhere and is the main 
term of the solution, whereas the second one is only im- 
portant close to the boundary, d F — z ~ I, when the 
exponents become of the same order. Each of the matrix 
functions hi (i = L, R) satisfies the equation 



A = 1 — 2i(uj — Ia z )r, 
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2hi(z) = Si(Xz) + - 



G{X{z - z'^h^dz 1 + - 
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G{X{z + z'))hi{z') + G(Xz')hi(z + z'))dz' , (13) 



where G{Xz) = G{Xz)e Xz / 1 and S L {Xz) = 
K(Xz) exp(Az). 

Far from the left boundary, z > I, we parameterize 
Hl{z) = A(z)l/Xz, 1/A = A -1 . Substituting it into 
Eq. (fl3|) and keeping the leading order in l/z, we obtain 
the equation for the diagonal matrix A, 
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2 - - In 2 A(z) = 1 
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A(z') 



dz' 



+ / A(z-z')G(Xz')dz', (14) 
XI J 

where the matrix £ = A J °° h,L(z)e~ 2Xz dz/l does not de- 



pend on z. The last term in Eq. ([14]) in the leading order 
in In -1 z is A(z)(j + ln(Az/7)) with 7 being Euler's con- 
stant. Subsequently, we obtain a differential equation for 
the function f* A(z')dz' j z' . The solution far from the 
boundaries reads 



h L {z)=A{z) 
A(z) = 



l e -\z/l 

~Xz~' 



6(a,X) 



(15) 



[2-(a/A)( 7 + log2Az/0]' 



where a constant 5(a, X) is of order one; at a = the 
exact solution gives 5(0, A) = 2. Numerical calculations 
show that 6 (a, A) is still close to 2 even for a = 1. 

Having solved the equation for we use it to deter- 
mine the matrix function Sr, according to 



Sr(z) = jj G(X(z + z')) 



h L (d F - z') - e- 2Xz '/ l h L (d F + z') 



dz' 



and find the solution for the function hn(z). 

Within the approximations used in the above analy- 
sis of the Eilenberger equation for the anomalous Green 
function /, the Josephson current density in the SFS 
structure can be represented as 
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Here, n(uj) is the Fermi distribution function. Substitut- 
ing the expressions for h(z) into Eq. (fTU|) and Eqs. ©, 
we find 
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where the matrix function Z depends on dp logarithmi- 
cally and for a = equals Z(Q, A, dp) = 1. Generally, in 




FIG. 1: Dependence of the factor \Z\ determined in Eq. fj 18p 
on the a — (tj — tv)/(tv + tj) for various values of It. 
From top to bottom: It = 0, 0.3, 1, 10. The results show that 
outside of the diffusive regime, /t<1 and tj > tv, Z is a 
smooth function of a of order of 1. 



the leading order in dp 1 it becomes 



Z =A^( 1+ A^ 



2 — — (1 + In 2) j A(dp), (18) 

where £ = X f Q °° liL(z)e~ 2Xz dz/l and £ = 
a~ 1 X 2 dp J °° hji(z)K(Xz)dz/l. For dp 3> I the quantity 
£ is constant, and A depends on it logarithmically. 
The quantity £ depends on dp in the same way as A. 
For At <C ||A|| ~ 1, the calculation of the frequency 
integral leads to an expression that is most conveniently 
represented as the sum over Matsubara frequencies 
uj n = 2-kT(ji + 1/2). This is equivalent to the replace- 
ment iuj — ■> oj n in the above expressions involving the 
matrix A. As the result, A becomes a diagonal matrix 
with two complex conjugate eigenvalues, so that in the 
Matsubara representation Z has the same property and 
can be written down as 

Z = \Z\ exp(ia z (pz) ■ 

Note that for uj n T < At < 1, \Z\ and (pz are two pa- 
rameters of a structure independent of the Matsubara 
frequency. For a = 0, one finds Z = 1. The dependence 
of \Z\ on the value of parameter a is plotted in Fig. [1] 
Although parameters \Z\ and ipz depend on the quanti- 
ties a = (tj — t\/)/(tv +tj), It, and dp/l, this fact does 
not qualitatively affect the results. Outside the regimes 
of It <C 1 and tj -C Ty (where our results are not appli- 
cable) Z is a smooth function of dp of the order 1 that 
does not contain any dependence on scales of order of 
£ rf ~ I or £ Q = vp/2I. 
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FIG. 2: The dependence of the critical current in SFS trilayer 
on ferromagnetic layer thickness. Here, we normalize current 
using a notional jo = 2-ireh'VFO 2 ''A, and show it for T = 0, 
At = 0.1 and several values of parameters a and It: a — 
(dashed lines) and a = 0.3 (solid lines). Dips in the value 
of j c indicate positions where it disappears and changes sign, 
thus resulting in a sequence of — n transitions. Neighbouring 
dips always correspond to the same value of It, demonstrating 
only weak dependence of the results on the parameter a. For 
comparison, we also show the decay of the Josephson prox- 
imity effect in an SNS structure heavily doped by magnetic 
scatterers (It = 0). 

Finally, we arrive at the expression for the critical cur- 
rent density [in j = j c sin((f>p — 4>r)\ which has the form 

e ~d F /l e -2uj Tl d F /v F 

cos {2I(1jp/vf + arctan(2ir) — (pz) 

In the limiting cases, the summation over Matsubara fre- 
quencies uj n can be calculated explicitly. For a ferro- 
magnetic layer with the thickness much greater than the 
coherence length in the superconductor, dp 3> vf/A, 
the sum equals 2T smh^^-nTdp /vp). In the oppo- 
site case of a thin layer, dp <C vp/A, one obtains 
(A/4) tanh(A/2T)). At zero temperature, the sum can 
be converted into integral which equals Ci(a) sin(a) + 
(ir/2 — Si(a)) cos(a), where a = 2Adp/vF, and the 
functions Si and Ci are sine and cosine integrals, re- 
spectively. For high temperature, dp 3> vf/T, only 



the lowest frequency is important, and the sum equals 
T(l + Tr^/A 2 )" 1 exp(-2nTdp/v F ). 

The critical current dependence on the ferromagnetic 
layer thickness described in Eq. (| 19|) for a weakly ferro- 
magnetic layer with dp > I is shown in Fig. [2] Even a 
when dilution of a ferromagnetic layer is such that ex- 
change energy in it is weak, It <C 1, oscillations of j c 
as a function of the layer thickness, with the period of 
£o = vf/2I, decay exponentially at the length scale of the 
mean free path, ^ = I - similarly to what happens in a 
disordered ferromagnetic layer with a strong exchange 10 
(It ^> 1). Our results for It 3> 1 coinscide with whose of 
Ref . fTol: For strong fields, the phase randomization of the 
order parameter is effective irrespective of the nature of 
scatterers. The dependence of the critical current on the 
thickness of the ferromagnetic layer in Eq. (fl"9)) resembles 
the experimentally observed suppression of the proximity 
effect, at the length scale comparable to the mean free 
path measured in the same material 2 -. Note that theories 
involving generation of the triplet order parameter due to 
nonuniform (spiral) magnetization in the ferromagnet 19 
end up with the opposite conclusion, predicting weaker 
decay of the order parameter. 

In conclusion, we developed theory of the proximity 
effect in a superconductor - weakly ferromagnetic GMR 
alloy - superconductor trilayers, which takes into account 
strong spin dependence of electron scattering of composi- 
tional disorder. The result, Eq. (fll?]) , describes 0-ir tran- 
sition for Josephson effect as a function of the thickness 
of the ferromagnetic layer dp: Oscillations occur with 
the period of £ Q = Vp /21 and exponential decay with the 
characteristic length ^ = I of the order of the mean free 
path, even in the regime when /t « 1. This result com- 
plements previous studies of the spin-singlet proximity 
effect in superconductor - ferromagnet hybrid structures 
performed for ballistic and diffusive systems with spin- 
independent scatterin g 6 i 10 i 13 as well as theories of the 
suppression of the order parameter oscillations caused 
by spin-active interfaces^ 

The authors acknowledge useful discussion with 
K.Efetov, the late A. Larkin, Yu. Nazarov, V. Ryazanov, 
and A. Volkov. This work was supported by EC Grant 
No. NMP2-CT2003-505587 and the Lancaster-EPSRC 
Portfolio Partnership No. EP/C51I743. 



1 T. Kontos et al, Phys. Rev. Lett. 86, 304 (2001). 

2 V. V. Ryazanov et al, Phys. Rev. Lett. 86, 2427 (2001). 

3 Y. Blum et al, Phys. Rev. Lett. 89, 187004 (2002). 

4 C. Bell et al, Appl. Phys. Lett. 84 , 1153 (2004); 
J. W. A. Robinson et al, |cond-mat/0606067| 

5 L. N. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, JETP 
Lett. 25 290 (1977). 

6 A. A. Golubov, M. Y. Kupriyanov, and E. Il'ichev, Rev. 
Mod. Phys. 76, 411 (2004). 

7 W. Guichard et al, Phys. Rev. Lett. 90, 167001 (2003). 



8 H. Sellier et al, Phys. Rev. B 68, 054531 (2003). 

9 T. Kontos et al, Phys. Rev. Lett. 89, 187004 (2002). 

F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 
B 64, 134506 (2001) 

1 A. I. Buzdin, L. N. Bulaevskii, and S. V. Panyukov, JETP 
Lett. 35, 172 (1982). 

2 A. A. Abrikosov and L. P. Gorkov, Sov. Phys. JETP 12, 
1243 (1961). 

3 A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

4 T. Yokoyama et al, Phys. Rev. B 71, 094506 (2005); 



■5 



L. Cretinon et al, Phys. Rev. B 72, 024511 (2005); 

M. Faure et al, Phys. Rev. B 73, 064505 (2006). 

S. Zhang, Appl. Phys. Lett. 61, 1855 (1992); S. Zhang and 

P. M. Levy, J. Appl. Phys. 73, 5315 (1993). 

T. Valet and A. Fert, Phys. Rev. B 48, 7099 (1993). 

G. Eilenberger, Z. Phys. B 214, 195 (1968); for an SFS 

junction see e.g. Ya. M. Blanter and F. W. J. Hekking, 



Phys. Rev. B 69, 024525 (2004). 

A. V. Zaitsev, Sov. Phys. JETP 59, 1015 (1984). 

See e.g., A. F. Volkov, A. Anishchanka, and K. B. Efetov, 

Phys. Rev. B 73, 104412 (2006). 

A. Cottet and W. Belzig, Phys. Rev. B 72, 180503 (2005). 



